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Abstract 

We discuss the stability and masses of topological solitons in QCD and strongly-interacting 
models of electroweak symmetry breaking with arbitrary combinations of two inequivalent 
Lagrangian terms of fourth order in the field spatial derivatives. We find stable solitons for 
only a restricted range of the ratio of these combinations, in agreement with previous results, 
and we calculate the corresponding soliton masses. In QCD, the experimental constraints on 
the fourth-order terms force the soliton to resemble the original Skyrmion solution. However, 
this is not necessarily the case in strongly-interacting models of electroweak symmetry break- 
ing, in which a non-Skyrmion-like soliton is also possible. This possibility will be constrained 
by future LHC measurements and dark matter experiments. Current upper bounds on the 
electroweak soliton mass range between 18 and 59 TeV, which would be reduced to 4.6 to 
8.1 TeV with the likely sensitivity of LHC data to the fourth-order electroweak Lagrangian 
parameters. 
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1 Introduction 



It is a generic feature of strongly-interacting theories that their vacua are not invariant 
under all their global symmetries, and that the resulting effective low-energy Lagrangians 
possess soliton solutions [l}[2]- The archetype is QCD, and many other four-dimensional 
examples have been studied theoretically [3). In the case of QCD, most studies have built 
upon the pioneering work of Skyrme (2] , and Skyrmion phenomenology has had a number of 
phenomenological successes (5j[6j . However, generalizations of the Skyrme model are possible 
even in QCD, and a fortiori in other strongly-interacting models that may exhibit dynamics 
rather different from QCD. 

Strongly-interacting models of electroweak symmetry breaking have long attracted con- 
siderable attention [71], with renewed intensity now that LHC experiments have discovered 
a Higgs-like particle |8l. However, although it has some of the characteristic properties pre- 
dicted for the Higgs boson in the Standard Model, the jury is still out, and many alternative 
scenarios described by low-energy effective chiral Lagrangians remain viable [9]. For example, 
there is active interest in the simplest possibility that the recently-discovered particle might 
be the pseudo-dilaton of some nearly-conformal strongly-interacting electroweak sector [To] . 
Discriminating between scenarios for electroweak symmetry breaking is a phenomenological 
priority, and the existence (or otherwise) of soliton solutions may be a valuable diagnostic 



tool for this task 11 



The masses and other properties of solitons in models described by effective low-energy 
chiral Lagrangians depend on the strengths of higher-order terms in their derivative expan- 
sions, and specifically of the coefficients of the fourth-order terms [2|. In the case of the 
minimal effective Lagrangian for SU(2) x SU(2) — > SU(2) that may be used to describe 
both QCD and electroweak symmetry breaking, there are two such parameters, as discussed 
in Section |2} One of these was considered originally by Skyrme p] , and has been the basis 
for most subsequent studies of Skyrmion phenomenology J5j. However, a priori the other 
term could also be present, and there have been some studies of generalized soliton solutions 



in the presence of this extra term 12 16 



Truncating the derivative expansion at fourth order is dictated primarily by practical 
reasons, since the contributions of higher-derivative terms to solitonic masses are not para- 



metrically suppressed 17 . Nevertheless, keeping only one specific fourth-order term, namely 
the Skyrme term introduced in Section [2j is phenomenologically satisfactory within 20 — 30%. 
Whether this accuracy is maintained with the inclusion of the non-Skyrme term can be tested 
in QCD, and we find that a somewhat larger range of classical masses is obtained in this case. 
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Since the calculated range of classical soliton masses overshoots the experimental nucleon 
mass in QCD, we infer that we can safely apply the truncated Lagrangian when looking for 
upper limits [j] 

The classical approximation is then followed by semiclassical quantization. Semiclassical 
rotations are suppressed in QCD in the 1/N C expansion 12]. So even if a given Lagrangian 
generates higher-order time derivatives, which is the case for the solitons with the non-Skyrme 
term present, they can be neglected within large N c approximation. A similar justification 
is applicable in many strongly-interacting models of electroweak symmetry breaking, but 
needs exploration on a case-by-case basis. Here we do not study issues beyond the classical 
approximation. 

As we review later in the context of our analysis, the existence or absence of stable soliton 
solutions hinges upon the ratio of the two fourth-order coefficients, and there is a generic 
range of this ratio where no stable solitons exist 12 -16 . Equally, there is a generic range, 
including the original Skyrmion as a special case, where solitons do exist and are stable, at 
least against spherically-symmetric decay. In the case of QCD, large-iV c arguments favour 
values of the fourth-order parameters within this stable range, close to the Skyrmion limit, 
and this possibility is also favoured by the available phenomenological estimates of these 



parameters in the effective chiral Lagrangian of QCD 18, 19 



However, currently we know very little about the possible magnitudes of the fourth-order 
coefficients in the electroweak case, and it is possible that their ratio is quite unlike the 
original Skyrme model, quite possibly in the range where no stable solitons exist. On the 
other hand, if the strongly-interacting electroweak sector is based on a theory with underlying 
constituents that bind to form 'electroweak baryons', one would expect solitons to exist and 
describe qualitatively their masses and other properties, even if they are rather different from 



the baryons and Skyrmions in QCD-like theories 11 



The purpose of this paper is to explore the non-Skyrmionic possibilities 12-16 , particu- 
larly in the electroweak case, discussing their existence, stability and masses in the classical 
approximation. This analysis is an essential ingredient in the exploration of the consistency of 
different strongly-interacting electroweak models with experiment. To this end, we assess the 



prospects for probing strongly-interacting electroweak models by confronting present 20 21 



and future limits [20] on (measurements of) fourth-order interaction parameters with limits 
on (measurements of) electroweak baryon masses. We also note that stable solitons would 
be present at some level in the Universe today as relics from the Big Bang and contribute to 

any case, one cannot expect that any phenomenological information on higher-order terms in the 
electroweak theory will be available in the very near future. 
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the dark matter [22] . One should check whether calculations of their abundance are compat- 
ible with cosmological and astrophysical estimates of the density of cold dark matter, and 



whether estimates of the rate for their elastic scattering on nuclei 23 are compatible with 
upper limits from direct searches for dark matter [24] . 

The structure of this paper is as follows. In Section [2] we discuss the mass and stability of 
an SU(2) soliton in the classical limit, including both possible fourth-order derivative terms. 
These calculations are used in Section [3] to set phenomenological bounds on the possible 
masses of QCD and electroweak solitons, taking into account the present and prospective 
experimental constraints on the fourth-order terms. Section [4] summarizes considerations 
concerning electroweak solitons as cold dark matter, and Section [5] contains some concluding 
remarks. 

2 Classical Mass and Stability of an SU(2) Soliton 

The effective chiral Lagrangian corresponding to some strongly-interacting sector with sponta- 
neously-broken SU(2) xSU(2) chiral symmetry is usually organized in terms of (even) powers 
of derivatives of the chiral field: 

U{x) = e X p(i^&) , (1) 

where the fields tt(x) correspond to the Nambu-Goldstone bosons (pions in the case of QCD, 
longitudinal polarization states of massive gauge bosons in the case of a strongly-interacting 
electroweak sector) and v is a typical symmetry breaking scale (v = F n = 93 MeV in 
QCD, v ~ 246 GeV in the case of electroweak theory). The derivative expansion is usually 
truncated at the fourth order, which could be reliable at energies below the characteristic 
strong-interaction scale, ~ GeV in QCD and ~ TeV or more in the electroweak theory: 

£ e ff = £2 + £4 , (2) 

where 

C 2 = ^Tr (dpUWrf) . (3) 

In the SU(2) case there are two independent invariants containing just four space-time deriva- 
tives (I2HI6I: 



£ 4 = 2s Tr [(R li R v )(R fl R v ) - (R^) 2 ]+2tTr [{R^R V ){R^R") + (R^f] , (4) 

where 

R, = d,UUK (5) 
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The parameters s and t can be in principle calculated from the underlying strongly-interacting 
theory by integrating out the constituent degrees of freedom (quarks and gluons in the case 
of QCD, and yet to be determined in an electroweak theory) or - in a phenomenological ap- 
proach - can be extracted from the data on the scattering of the Nambu-Goldstone bosons 
or massive gauge bosons (20) In QCD, as we discuss in more detail later, the large- N c 
expansion and pion data indicate that \t\ <C \s\, so that the effective Lagrangian contains 
(in a first approximation) only the kinetic term ^ and the first term in Q, as in Skyrme's 
original work. For this reason, we refer to this as the Skyrme term, and refer to other term 
in Q, that with coefficient t, as the non-Skyrme term. 

Over fifty years ago already, Skyrme |2j observed that (J5]) with t = possesses solitonic 
solutions that, due to the fact that 7r 3 (SU(2)) = Z, carry an integer-valued topological 
quantum number 

B = i / d3xel ' kTl [(U^d^^d^^dkU)] . (6) 

that can be interpreted as baryon number (see also (l)). 

It is straightforward to calculate the classical contribution to the mass of the baryon in 
terms of the parameters v , s and t. In principle, all terms in the derivative expansion make 
significant contributions to the soliton mass, whereas we have no information on the possible 
magnitudes of higher-order terms. On the other hand, we know that classical calculations 
in QCD keeping only the fourth-order term Skyrme term, i.e., setting t = 0, are accurate 
to ~ 20 — 30% [5lp3|. In the case with non-zero non-Skyrme term a somewhat larger masses 



are allowed, cf. our calculation below for QCD point A in Sect |3.1| We therefore hope 
that classical calculations in the electroweak case, keeping the general form of fourth-order 
coupling, sample the range of possible masses with similar accuracy. 

Refinements to this classical calculation include calculation of the Casimir energy, inclu- 
sion of the current-quark mass terms and semiclassical quantization of the rotational modes, 
and provide better agreement with the observed baryon masses. However, we do not go into 
them here, for three additional reasons, in addition to our ignorance of possible higher-order 
terms. One is that we are primarily interested in bounds on the masses of the possible elec- 
troweak baryons, rather than in details of their spectrum. A second reason is that, in any 
case, we lack the information about the underlying strongly-interacting theory that would 
be needed to calculate the non-classical corrections in this case. The third reason is that 
in reality, even in QCD, the parameter t is not necessarily negligible. Instead, our primary 



2 For rigorous lower bounds from Lorentz invariance, analyticity, unitarity and crossing, see 



21 
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interest here will be the classical soliton mass in the presence of such a non-negligible non- 
Skyrme term, and how the possibility that t ^ affects the (approximate) mass bounds 
provided by the Skyrme calculation with t — discussed in [II]. 

Classical soliton solutions of the chiral field equations are usually found within the 
spherically-symmetric 'hedgehog' Ansatz for a static field configuration: 

U{f) = exp ^-^P(r)J , (7) 

where the profile function P(r) is a solution to the Euler-Lagrange equation of motion, which 
for the Lagrangian (^-Q takes the following form: 

P"W = -fg (8) 

P(0) = TT, P(oo) =0. (9) 



with boundary conditions 



Here 

F(p) = p 2 (l - 96 tP' 2 ) + 64 s sin 2 P . 



= -64 tpP' 3 + 32 s sin 2PP' 2 + 2pP' 

in 1 
P 



sin 2 P N 

sin 2P ( 1 + 32 (s - t) ) , (10) 



and the dimensionless radial variable p is defined as 

p = rv , (11) 

where v = f n , ~ 246 GeV in the QCD and electroweak cases, respectively 

Following the pioneering paper |2j , classical solutions of the field equations were studied 
numerically long ago in [12] and stability regions in the (s, t) parameter space were analyzed 



in 



13 -16 . We have repeated their analyses for the purpose of the present work with the 



following results. 

The total classical soliton mass 



M sol = Ativ (M 2 + M 4 ) (12) 



where 



M 2 

2 



oo 

^ dp (p 2 P' 2 + 2 sin 2 P 



o 

oo oo 



M 4 = 16s J dp [2P 12 sin 2 P + ^f-j -Stjdr (^p 2 P u + 2^f- ) , (13) 
o 



5 



where the contribution corresponds to £2 and is clearly always positive, and the contri- 
bution M 4 corresponds to £4. By changing variable once more: p — > pj \f\s\ or p — > pj ^/|t[, 
one can show that the solutions depend only on the ratio t/s. By rescaling: p — > Xp, one can 
show that solutions exist only when M4 is positive. Inspecting directly M 4 in (13), we see 
that in the fourth quadrant of the (s,t) plane (s > 0, t < 0) M4 is always positive, so that 
solutions of (|sj) always exist, whereas in the second quadrant (s < 0, t > 0) M 4 is always 
negative, so there are no solutions of d8j). In the first and third quadrants positivity bounds 



can be derived 12, 13 . It turns out that in the first quadrant M4 is negative for t/s > 2, 
so the soliton is unstable, and in the third quadrant M4 stays positive for t/s > 2, so the 
soliton is stable. This is illustrated in Fig. [TJ where the regions of positive and negative 
M 4 are displayed by green shading (diagonal squares) and red shading (vertical squares), 
respectively. 



I instability 




non-Skyrmc branch 



Figure 1: The (s,t) parameter plane, indicating with green shading (diagonal squares) the 
range of parameters satisfying the positivity bound M 4 > and hence admitting a solution of 
the differential equation and with red shading (vertical squares) the range where solutions 
do not exist due to the negativity of the fourth order term: M 4 < 0. In the range t/s < 0.29 
within the first quadrant, indicated in blue (vertical bars), solutions with positive M 4 exist, 
but no positivity bound can be derived. 

In the remaining parts of quadrants I and no positivity bounds can be derived and 
the positivity of M4 has to be checked numerically. It turns out that in the third quadrant 
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M so \ > and — > + when t/s — > 2 + , whereas no solution exists for t/s < 2. The situation 
is different in the first quadrant, where solutions with positive M4 exist for small t/s below 
the line t/s ~ 0.29, where the function F(p) vanishes and equation ^ cannot be solved. 
Interestingly, F(p) is equal to the second variation of M so i, so the positivity of F(p) is a 
necessary condition for M so \ to be the minimum with respect to the variations of P(r). 
Below the line t/s ~ 0.29 the solution is classically stable, and above this line it is classically 
unstable. 



However, it has been shown in 12 that there is no lower limit on the soliton mass 
anywhere in the first quadrant. This means that, when t/s < 0.29, the classical solution 
discussed above, though stable against local spherically-symmetric perturbations, can only be 
metastable at best, and might be unstable against non-spherically-symmetric perturbations. 
This region is shaded blue (with vertical bars) in Fig. [I] 

As already remarked, soliton solutions depend only on the ratio t/s (or s/t). We analyze 
solutions starting from the Skyrme-like case with t = and the non-Skyrme-like case with 
s = 0, and have divided the allowed parameter space into two branches on either side of the 
line t/s = —1, as also depicted in Fig. [TJ In order to solve equation ^ within the Skyrme 
branch where s 7^ 0, we have rescaled p — > p/y/s and evaluated the soliton mass in units of 
47rv\/s, whereas within the non- Skyrme branch where t 7^ we have rescaled p — > p/\f— t 
and evaluated the soliton mass in units of Atcv t. The results are presented in Fig. [2] as 
functions of e = t/s for the Skyrme branch and e = s/t for the non-Skyrme branch. 

Finally, in Fig. [3] we present a contour plot of the classical soliton mass in units of 4irv 
over the allowed range of parameter space. We see that the contours of constant M so \ are 
almost parallel to the line t/s = 2 where the soliton mass vanishes. This feature is very 
helpful for finding the maximum of the generalized soliton mass, as we discuss in the next 
Section. 

3 Phenomenological Estimates of Soliton Masses 

We now discuss constraints on the conventional baryon mass using constraints on the param- 
eters of the low-energy effective chiral Lagrangian for QCD, and constraints on the masses 
of possible 'electroweak baryons' given by the constraints on higher-order electroweak La- 
grangian parameters from electroweak theory, current phenomenology and the potential 



sensitivity from the LHC 11,20 



Phenomenological constraints on higher-order Lagrangian parameters are often given in 
terms of the coefficients 0^5 that are related to the parameters s and t discussed in Section [2] 
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Figure 2: The soliton mass in units ofAnv^s for the Skyrme branch where e = t/s (solid 
blue line), and for the the non- Skyrme branch in units of 47rvy/^i where e = s/t (dashed red 
line). The Skyrme branch ends at t/s ~ 0.29 (see Fig. [7]). 

in the following way [TT] : 

s = t = (14) 

Although limits on electroweak baryon masses are our principal interest, we first discuss the 
limits that can be derived for ordinary baryons in QCD, comparing the non-Skyrmion case 
with t 7^ with the conventional Skyrmion case t — 0. 

3.1 QCD baryons 

Higher-order coefficients in the low-energy effective Lagrangian of QCD have been the sub- 
ject of many studies. Here we use the following ranges extracted from low-energy strong- 
interaction data: [l~l~|[l8|[l~9] : 

llxl(r 4 < a 4 <17xl(r 4 , 

14 x 10" 4 < a 4 - « 5 < 40 x 10" 4 . (15) 
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Figure 3: Contour plot of the mass of the generalized Skyrmion in units ofAwv. 

We can compare these values with the predictions of the large- N c approximation within the 
framework of chiral SU(3) x SU(3) SU(3) pi: 



(16) 



a 4 = 18 x 10~ 4 , 
a 5 = -16xl0 -4 . 



The bounds ( fl5"| are displayed in Fig. |4| where the large- N c prediction is also indicated, as 
a red spot f\ 

As already discussed in Section [2] the contours of constant M sol are almost parallel to the 
line t — 2s. Therefore the maximal mass corresponds to the right-most corner of the region 
allowed by the QCD bounds, called Point A below. For completeness, we also include two 
pure Skyrme points (B and C) and the point corresponding to minimal mass, located right 
below the stability line t/s = 0.29 (Point D): 

• A : (s,t) = (10, -4.5) x 10~ 4 (maximal mass point overall): M A ~ 1354 MeV; 

• B : (s, t) = (8.5,0.0) x 10 -4 (maximal mass Skyrme point): M B ~ 1118 MeV; 



3 We represent the constraints (151, ( fl9[ ) , (20) as parallelograms in the (s,t) plane, whereas they should 
be ellipses. However, the correlations between the errors in s and t are not available, so these ellipses cannot 
be drawn accurately. 
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• C : (s,t) 

• D : (s,t) 



= (5.5,0.0) x 10 4 (minimal mass Skyrme point): M c ^ 900 MeV. 
= (4.3, 1.2) x 10" 4 (minimal mass point overall): M D ~ 728 MeV. 




Figure 4: Comparison of the parameter ranges allowed for the low-energy chiral perturbation 



theory coefficients in QCD '18, 191, eq. (15), with the range of parameters where stable 



baryonic solitons exist. The parameters a^, a^, s and t are in units of 10 4 . The red 



spot corresponds to the large-N c values given in eq. (16). Point A has the maximal mass 
overall, Points B and C bracket the range of parameters allowed in the pure Skyrme limit: 
t = 0, and Point D lies on the stability boundary t/s = 0.29 and has the minimal mass 
overall. 

We see that the range of classical masses allowed in the Skyrme case t = includes the 
physical value of the nucleon mass. In terms of the conventional representation 

' (17) 



32e 2 ' 

this range corresponds to 

6 < e < 7.5. (18) 

We also note that the range of possible QCD baryon masses is extended significantly in 
the presence of a non-Skyrme term, by (9(200) GeV in either direction. However, detailed 
QCD baryon phenomenology including the evaluation of corrections to the classical soliton 
mass due to the Casimir energy, inclusion of the current-quark mass terms and semiclassical 
quantization of the rotational modes lies beyond the scope of this paper. 
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3.2 Current bounds on elect roweak baryon masses 

The higher-order coefficients in an effective electroweak Lagrangian are currently poorly 
constrained, namely by the following bounds |20|: 

-3.5 x 1CT 1 < a 4 < 0.6 x 10" 1 , 

-8.7 x 10" 1 < ct 5 < 1.5 x HT 1 . (19) 

These constraints are superimposed on the (s, t) plane in Fig. [5j As in the QCD case, we 
superpose various points that illustrate the range of possible electroweak baryon masses 
currently allowed, as obtained using the values of the soliton mass calculated in Section [2] 
and displayed in Fig. [2} We note that there is no minimal electroweak baryon mass, since the 



constraints (19) include the possibility that M so \ = for t = 2s < 0. The points displayed 
include the maximal mass point denoted by A, the maximal-mass Skyrme case (Point B), 
and the maximal-mass case with vanishing Skyrme term s = (Point C): 



• A : 




(0.23, 


—0.20) (maximal mass point overall): 


M A 


~ 59 TeV; 


• B : 


(M) = 


(0.03, 


0.0 ) (maximal mass Skyrme point): 


M B 


~ 18 TeV; 


• C : 


M) = 


(0.0, - 


-0.175) (maximal mass with s = 0): 


M c 


~ 31 TeV. 



This analysis indicates that the upper bound on the possible electroweak baryon mass 
over the range of parameter space currently allowed is approximately 60 TeV This upper 
limit is somewhat higher than in |TTJ, because we have extended the analysis to the case 
of a non-zero non-Skyrme term. If the conjectural strongly-interacting electroweak theory 
were to resemble QCD, i.e., the parameter t were small, the upper limit on the possible 
electroweak baryon mass limit would be significantly reduced to ~ 18 TeV. 

3.3 Prospective LHC bounds on electroweak baryon masses 

A study has been made of the prospective LHC sensitivity to the higher-order electroweak 
Lagrangian terms. It was estimated that, in the absence of a signal, the LHC could yield 



the following allowed ranges 20 



-7.7 x 10 -3 < a 4 < 15 x 10~ 3 



-12 x KT 3 < a 5 < 10 x 10~ 3 . (20) 



'Assuming that the model-dependent corrections to the classical mass calculation are not very large. 



11 




Figure 5: Comparison of the parameter ranges currently allowed t 20^ for the effective low- 
energy electroweak Lagrangian with the range of parameters where stable electroweak baryonic 
solitons exist. The parameters a^, as, s and t are in units of 10 _1 . Point A has the maximal 
mass overall, Point B has the maximal mass allowed in the pure Skyrme limit: t = 0, and 
Point C has the maximal mass allowed in the limit where the Skyrme term vanishes: s = 0. 



These ranges are superimposed on the (s, t) plane in Fig. |6j As in the previous examples, we 
list below the classical masses calculated for three illustrative points that are also displayed 
in Fig. [6] 

• A : (s,t) = (6.75,0.75) x 1(T 3 (maximal mass point): M A ~ 8.1 TeV; 

• B : (s, t) = (6.0, 0.0 ) x 10~ 3 (Skyrme point): M B ~ 7.9 TeV; 

• C:(s,t) = (0.0, -3.85) x 10~ 3 (non-Skyrme point): M c ~ 4.6 TeV. 

The upper limit on the possible electroweak baryon mass in the event that the LHC does 
not find non-zero higher-order effective Lagrangian parameters is ~ 8 TeV, very similar to 
the bound quoted in [11] . This reflects the fact that the point with maximal mass overall is 
quite close to the Skyrme limit t = 0. In this case, unlike the previous examples, the upper 
limit on the soliton mass is not relaxed by allowing a non-Skyrme term with t < 0. 
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Figure 6: Comparison of the estimated LHC sensitivity range f2uy for higher-order elec- 
troweak Lagrangian parameters with the the range of parameters where stable electroweak 
baryonic solitons exist. The parameters a^, «5, s and t are in units o/10 -3 . Point A has 
the maximal mass overall, Point B has the maximal mass allowed in the pure Skyrme limit: 
t = 0, and Point C has the maximal mass allowed in the limit where the Skyrme term 
vanishes: s — 0. 

4 Electroweak Baryons as Dark Matter? 

If they exist, electroweak baryons should be present in the Universe today, and could provide 



cold dark matter 22 . Their possible relic density today depends on the presence and mag- 
nitude of a primordial electroweak baryon asymmetry. If this was small, and the electroweak 
baryon density was completely equilibriated in the early Universe, the present abundance of 
electroweak baryons would be insufficient to provide the present density of cold dark mat- 
ter. In this case, one or more other sources of cold dark matter would be required, and the 
constraints on dark matter scattering that we now discuss would be irrelevant. 

Assuming that they make up the bulk of the cold dark matter density, it is clear that 
electroweak baryons should have no electric charge. Moreover, it was pointed out in [25] 
that they could not be fermions, as these would have an unacceptably large cross section 
for spin-independent scattering through magnetic moment couplings. Some scenarios for 
the effective low-energy electroweak Lagrangian in which the lightest electroweak baryon is 



a neutral boson have been enumerated in 26 , in an analysis based on the Wess-Zumino 
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Witten term and the standard relation Q cm = I3 + Y. The neutral boson scenarios found 
do not include SU(3) x SU(3) — >■ SU(3) (which could yield a neutral fermion if iV c = 3 or a 
bosonic baryon if N c is even), but do include SU(iV) x SU(iV) SU(iV) with N > 3 and 
N c even, and SO(iV) x SO(iV) -> SO(iV). However, SU(iV) ->■ SO(iV) would yield a boson 
with charge N c , and the other coset structures considered in |26| have trivial homotopy and 
hence no electroweak baryons. 

It should be noted that dark matter scattering is a non-trivial constraint, even if the 
lightest electroweak baryon is a neutral boson. The scattering cross section for this case was 



estimated in 23 , with results indicating that models containing a pseudo-dilaton identified 
with the recently-discovered boson X with mass ~ 125 GeV would be consistent with the 
(updated) XENON100 limit 24 if the pseudo-dilaton couplings were scaled by > 1 TeV Q 



and the electroweak baryon weighed < 1 TeV, as seen in Fig. 4 of [23]. The XENON100 
Collaboration does not report results for larger dark matter particle masses, but naive ex- 
trapolation of the theoretical calculations and experimental sensitivity would suggest some 
tension for masses above 1 TeV. 

Our results suggest that scenarios which, according to the analysis of [26], might yield ei- 
ther a charged bosonic baryon or a fermionic baryon should not necessarily be abandoned. By 
the same token, one might not need to worry about the prospective tension with XENON100. 
This is because the parameters of the fourth-order effective Lagrangian might be in the range 
where no stable solitonic baryon exists, i.e., t > or t < and t/s > 2. That said, it could 
be that stable electroweak baryons nevertheless exist for such parameter choices, but cannot 
be found using the solitonic Ansatz discussed here. 

5 Concluding Remarks 

In this paper we have revisited the existence, stability and possible masses of classical soli- 
tonic solutions to the effective low-energy Lagrangians of QCD and a possible strongly- 
interacting electroweak sector, with particular attention to the presence and consequences 
of a possible non-Skyrme quartic term as in (|4]) 12 14 16 . We have revisited the stability 



constraints in the (s, t) plane, see Fig. [TJ and given general results for the classical soliton 
mass, see, e.g., Fig. [3j 

In the QCD case, we have found that current phenomenological constraints on the quartic 



s This is in the ball-park calculated in holographic dilaton models, and is reported to be consistent with 
the available data on X(\2h) [27], once the contributions of the strongly-interacting electroweak sector the 
Xgg and X77 couplings are taken into account. 
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chiral Lagrangian terms 18,19 allow a range of classical baryon masses that is somewhat 
broader than in the pure Skyrme case with t = 0. In the case of a possible strongly-interacting 



electroweak sector, current data 20 allow a somewhat wider range of t/s than in QCD, and 
the expansion in the possible range of electroweak baryon masses is proportionally larger 
than in the pure Skyrme case. On the other hand, the prospective LHC sensitivity [20] 
to quartic Lagrangian parameters would not allow masses substantially larger than in the 
Skyrme limit. 

It would be interesting to establish a 'no-lose' theorem that the LHC will either find 
non-zero higher-order Lagrangian parameters or discover electroweak baryons. However, 
this would require greater LHC sensitivity than indicated by the prospective range (20). To 
our knowledge, there is no accurate estimate of the range of electroweak baryon masses that 
the LHC could detect, but it is surely less than the ranges quoted above on the basis of a 
theoretical estimate of the possible LHC sensitivity to higher-order Lagrangian terms. On 
the other hand, the prospective LHC sensitivity should be revisited by a full experimental 
simulation including the possibility of high-luminosity LHC running. 

As already commented, in order to be a valid dark matter candidate, a stable electroweak 
solitonic baryon should be a neutral boson. Even in this case, Fig. 4 of ^3] indicates that such 
a model is quite constrained, at least in models in which the recently-discovered new boson 



is interpreted as a pseudo-dilaton. The estimate given in 23 of the spin-independent cold 
dark matter scattering cross section rises with the Skyrmion mass, so that experiments such 
as XENON100 [24] may be more sensitive to models with heavier electroweak baryons. Full 
exploration of this issue lies beyond the scope of this paper, but it is clear that dark matter 
scattering experiments are potentially interesting probes of strongly-interacting electroweak 
models with Skyrmion solutions. 
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